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performed for aeroservoelasticity (2.1.7) tasks under the Technology 
Integration (2.1) effort. Boeing personnel performed the analysis work reported in 
this document. 
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Summary 


The HSCT Flight Controls Group is developing a longitudinal control law, known 
as Gamma-dot / V, for the NASA HSR program. Currently, this control law is 
based on a quasi-steady aeroelastic (QSAE) model of the vehicle. This control 
law was implemented into the p-k flutter analysis process for closed loop 
aeroservoelastic analysis. The available flexible models, developed for the TCA 
aeroelastic analysis, were used to assess the effect of control laws on flutter at 
several different Mach numbers and mass conditions. 

Significant structures and flight control system interaction was observed during 
the initial assessment. Figures 1 and 2 present a summary of the effect of total 
closed loop gain and phase on flutter mechanisms, based on ideal sensors and 
real sensors, for Mach 0.95 and mass M02 condition. Control laws based on 
ideal sensors gave rise to increased coupling between the rigid body short period 
mode and the first symmetric elastic mode. This reduced the stability margins 
for the first elastic mode and does not meet the required 6 dB gain margin 
requirement. The effect of “real” sensors significantly increased the structures 
and control system interactions. This caused the elastic, modes to be highly 
unstable throughout most of the flight envelope. 

State-space models were developed for several conditions and then MATLAB 
program was used for the aeroservoelastic stability analysis. These results 
provided an independent verification of the p-k flutter analysis findings. Good 
overall agreement was observed between the p-k flutter analysis and state- 
space model results for both damping and frequency comparisons. These results 
are also included in this document. 

The following conclusions were made during this study: 

1) The traditional QSAE control design strategies are not applicable; control 
laws need to be developed based on elastic models. 

2) The location of sensors impacts various elastic modes differently. 

3) A MIMO design strategy may be required. This would add complexity to the 
design architecture and raises the challenge for robust design. 

4) The aeroservoelastic analysis using p-k and state-space formulations 
provided almost identical results. 

5) The p-k based flutter analysis and state-space formulation analysis methods 
should be used concurrently for design and verification. 


1 6 1 . ,acct;KSfi8042;ASE/M95/M02/t)595l5a6bgv2.esb 
C] : /acc t /k s n8042 / ASE / M95 / M02 / b595 1 5a6b p v 1 . esb 

[D] : / ac c t / k sn8042 / ASE/M95/M02/ b595t 5a6cpv6 .esb 

[E] : / ac c t / k s n 804 2 / ASE / P LOTS / GV / g vma c h 95 .esb 


cp pp-cl bQie.' 


C*LC k.S.NAGARAJA 26Mar98 


CONTROL SYSTEM EFFECTS 
STRENGTH+FLUTTER SIZED AIRPLANE 
PITS MODEL TCAY 

JB£FJEJAf£r 





















\ 1 : / acc t /k sn8 042/ ASE /PLOTS/ GV/gvmach95 . esb 

[B] : /acc t / ksn8042/ASE/ M95 / b 595 1 5a2gv 1 . esb 

[C] : / ac c 1 / k sn8042/ ASE /M95/ b595 15a2pv1.esb 



cue K.S.NAGARAJA 26Ma r 98 


CONTROL SYSTEM EFFECTS 
STRENGTH+FLUTTER SIZED AIRPLANE 
DITS MODEL TCAY 
























Introduction 


This report presents work performed by The Boeing Company to satisfy the 
deliverable “Evaluation of aeroservoelastic Effects on Symmetric Flutter” for 
Subtask 7 of Reference 1 . The objective of this report is to incorporate the 
improved methods for studying the effects of a closed loop control system on the 
aeroservoelastic behavior of the airplane planned under NASA HSR technical 
Integration Task 20 work. Also, a preliminary evaluation of the existing pitch 
control laws on symmetric flutter of the TCA configuration was addressed. The 
goal is to develop an improved modeling methodology and perform design 
studies that account for the aero-structures-systems interaction effects. 

Approach 

Aeroelastic Models 

Aeroelastic models developed for the TCA aeroelastic analysis, discussed in 
Reference 2, were considered for this study. Appendix A summarizes the 
analytical approach used for the control law assessment study. It presents the 
analysis process, Elfini FEM based dynamics models and aerodynamic mesh for 
the subsonic and supersonic linear aerodynamics. The dynamic math model 
utilized 60 flexible modes, 3 symmetric rigid body modes, and two assumed 
degrees of freedom for the stabilizer and elevator rotation modes required for 
control system closed loop representation. The mass cases considered for the 
study and the required flutter clearance speed envelope (Vd / Md and 1 .15 Vd / 
Md) are also provided in the Appendix A. The elastic models and the linear 
unsteady aerodynamic theories used for the aeroelastic analysis are discussed 
in Reference 2. The strength and flutter sized stiffness designed model 
(identified as TCAY) was used to assess the effect of control laws on flutter at 
several Mach numbers and mass conditions. 

Control Laws 

The HSCT Flight Controls Group is developing a longitudinal control law, known 
as Gamma-dot / V, for the NASA HSR program. Currently, this control law is 
based on a quasi-steady aeroelastic (QSAE) model of the vehicle. The 
Controller has an inner loop (SAS system), an outer loop (control path) and a 
feed forward direct control path. The command inputs to the stabilizer and 
elevator actuators are functions of flight path angle, angle of attack, pitch 
attitude, and pitch rate measured at the IMU location (BS1211). The various 
gains in the control law are schedule to be functions of dynamic pressure. The 
control law was transformed into inertial reference axis before implementation 
into the p-k flutter analysis. Also, the analysis utilized a simplified third order 
actuator model at the present time. The Apex based p-k flutter equations of 

4 






motion from Reference 3, simplified control law, sensor equations, and the 
actuator model used for the closed loop flutter analysis are summarized in 
Appendix A. The details of the longitudinal control law are presented in 
Appendix B. 

Results 

Flutter Analysis based on Apex p-k Solution Process: 

Matched point flutter analyses were performed, for Mach numbers .95 and 2.6 
and mass conditions M02, MT1, and MT4, to get a better assessment of th'e 
effect of longitudinal control law on flutter stability. Significant interaction 
between the structural modes and control system was observed in the low 
frequency range of up to 5 Hz. Initial assessment was focused on this frequency 
range and the results are summarized below. 

Mach 0.95, Mass case M02: 

Open loop analysis : Figure 3 shows that all the aeroelastic modes are stable in 
the speed range of interest. 

Closed loop analysis based on ideal sensors (Mean axis feedback) : This 
assumes that the sensors at the IMU location respond to rigid body motion of the 
airplane and do not respond to the elastic motion of the structure during the 
oscillations. Figure 4 shows the flutter solution that includes the effect of the 
control system for the nominal gain of the control law. The results show that 
there is increased coupling between the rigid body short period mode and the 
first elastic structural mode (1.45 Hz). This increased coupling contributed to the 
marginally unstable behavior of the first elastic mode. This does not meet the Vd 
flutter clearance requirement for zero structural damping. Figure 5 shows the 
flutter mode shape for this flutter mechanism. This exhibits significant elastic 
motion of the outboard nacelle and the wing combined with aft body vertical 
bending motion. Figure 6 shows the control system gain and phase variation 
effects on flutter speed, with a structural damping value of 0.03, and the gain and 
phase margins for the critical modes. These results show that the 
aeroservoelastic system does not meet the required 6 dB gain margin for the 
1 .45 Hz elastic flutter mode. 

Closed loop analysis based on real sensors (Full structural feed back analysis) : 
This analysis accounts for the real structural deformation feed back during the 
oscillations. Figure 7 shows the flutter solution for this situation and the control 
laws caused the 1 .7 Hz flutter mode to become highly unstable in the speed 
range of interest. Figure 8 shows the flutter mode shape for this flutter 
mechanism and again exhibits significant elastic motion of the outboard nacelle 
and the wing combined with the whole body flexing of the airplane. Figure 9 
shows the effect of system gain and phase variation effect on flutter speed. The 
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flutter mechanism is highly unstable and a significant level of gain attenuation is 
required to stabilize the flutter mode at this frequency either by reducing the 
system gain or by using a deep notch filter. This would result in a significant 
penalty on airplane performance. The IMU pitch rate feed back of the stability 
augmentation system was identified as the main cause for the increased 
structural / system interaction. An effort was made to study the effect of sensor 
location on the aeroservoelastic stability. Figures 10-13 show the effect on 
flutter speed based on sensors located near C.G. of the airplane. The results 
based on ideal sensors still exhibit low unstable damping for the 1 .45 Hz mode. 
Also, the results based on real sensors (full structural feed back) show a different 
flutter mechanism (3.6 Hz) that is again highly unstable and does not meet 6 dB 
gain and 60 degree phase margin flutter requirements. 

Mach 2.60, Mass case M02: 

Open loop analysis : Figure 14 shows that all the aeroelastic modes are stable in 
the speed range of interest. 

Closed loop analysis based on ideal sensors ; Figure 15 shows a marginally 
unstable 1 .6 Hz mode. This does not meet the Vd flutter clearance requirement 
for zero structural damping. However, Figure 16 show that the flutter mechanism 
exhibits adequate gain and phase margin with structural damping of 0.03. 

Closed loop analysis based on real sensors : Figures 17 and 18 show flutter 
characteristics similar to Mach 0.95 solution discussed above and the results do 
not meet flutter clearance requirements. Again, the effect of sensor location 
(near C.G.) on flutter characteristics was examined and the results shown in 
Figures 19 through 22 demonstrate flutter characteristics very similar to the 
Mach 0.95 solution. 

Mach 0.95,2.60 for Mass cases MT1 and MT4: 

Figures 23 to 42 show the flutter analysis results for Mach numbers 0.95 and 
2.60 and Mass cases MT1 and MT4. These results include the open loop 
solution, closed loop results including gain and phase variations based on mean 
axis analysis (ideal sensors), and full structural feed backs (real sensors). The 
results show the flutter characteristics to be very similar to the cases discussed 
above. 

Stability Analysis based on State-Space Formulation Process: 

The above aero-structures-systems interaction effects were independently 
validated by developing state-space models for several Mach / Mass conditions. 
The methodology used for this process is described in Appendix C. The analysis 
utilized 60 flexible modes, two rigid body modes, and two assumed stabilizer 
rotation and elevator rotation modes. Also, a sensitivity analysis, for number of 
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aerodynamic lag (terms varying from 4 to 8, was performed. For comparison of 
results with p-k flutter analysis results it was decided to use 8 lag terms in the 
state-space formulation process. The stability analyses were performed for each 
of the models at various flight speed conditions using MATLAB program. 

Figures 43 through 46 present the open loop root locus match point stability 
solution results for M02 and MT1 mass cases at Mach numbers 0.95 and 2.60. 
These results show that the aeroelastic modes are stable at all speeds 
considered in the 1.15 Vd envelop thereby validating the p-k open loop flutter 
solutions discussed earlier. 

Figures 47 and 48 present the mean axis root locus match point stability solution 
results (ideal sensors) for M02 and MT1 mass cases at Mach numbers 0.95. 

The results show that the low frequency elastic mode is marginally unstable for 
nominal system gain predicting characteristics similar to the p-k flutter analysis 
results. 

Figures 49 through 52 presents the root locus match point stability solution 
results based on full structural feed backs (real sensors) for the above mass 
case / Mach number combinations. The low frequency flutter mechanism exhibit 
high instability as predicted by the p-k flutter analysis process. In addition, there 
are additional high frequency mechanisms (8.0 to 9.5 Hz) that are highly 
unstable. Further work is needed to understand these high frequency flutter 
mechanisms. 

Comparison of P-K Flutter and State-Space Model Analyses Results : 

Mach 0.95 and 2.60 with mass case M02 conditions were considered for 
quantitative comparison of the stability analyses results of the two processes 
considered above. This was done by transforming the state-space model results 
into the frequency / damping format for comparison with p-k frequency / damping 
results. Figures 53 and 54 present the match point open loop comparison 
between the two sets for Mach numbers 0.95 and 2.60 respectively for M02 
mass case. The results show good comparison of frequency between the two 
processes and acceptable level of damping comparison. 

Figure 55 show similar results for the mean axis closed loop analysis for the 
mass / Mach combination case M02 / 0.95. Both the processes predicted similar 
level of low damping instability for the 1 .45 Hz elastic mode. However, there is a 
significant difference in the short period mode frequency at higher speeds 
between the two methods. Further investigation would be done to understand 
this difference. 

Figures 56 and 57 show results of comparison based on real sensors (full 
structural feed backs) analysis. The results show good comparison of 
frequencies between the two processes and acceptable level of comparison for 
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stable modal damping values. Also, the highly unstable damping value 
comparison between the two methods, for the low frequency mode 1 .45 Hz, is 
considered to be good. 


Conclusions 

This aeroservoelastic investigation has demonstrated the following: 

1) The QSAE based longitudinal gamma-dotA/ control law drastically violates 
the flutter stability requirements through most of the flight envelope. 

2) Traditional QSAE control design strategies are not applicable. 

3) Control laws need to be developed based on flexible dynamic models. 

4) The location of sensors impact various elastic modes differently. 

5) The aeroservoelastic analysis using p-k and state-space formulations 
provided almost identical results. 

6) The p-k based flutter analysis method and state-space formulation analysis 
method should be used concurrently for design. 

7) Early aeroservoelastic analysis for ASE wind tunnel model is required for 
defining control laws, and may impact wind tunnel model design. 


Recommendations 


1 ) Assess the effect of flexibility on other control laws such as flare, ride 
qualities and lateral/directional control. 

2) Incorporate Long Beach P-transform state-space modeling approach into the 
Flight Controls / Structures tool suite. 
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This page will be updated at a later time to include p-k flutter solution for 
M=2.60, Mass Condition: MT1, Open Loop Solution 
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This page will be updated at a later time to include p-k flutter solution for 
M=2.60, Mass Condition: MT4, Open Loop Solution 

Figure 38 
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M=2.60, Mass Condition: MT4, Gain & Phase Variation (Mean Axis) Solution 
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APPENDIX A 


Analytical Approach for Closed Loop 
Flutter Analysis 


Task 20: Sub Task 7 - Aeroservoelastic Design Studies 


High Speed Civil Transport 

ANALYSIS PROCESS 


• STRUCTURAL DYNAMIC MODEL : DITS TCAY MODEL WITH STIFFENED 
ACTUATOR SPRINGS FOR STABILIZER AND ELEVATORS 

• PITCH CONTROL LAWS (GAMMA DOT - V CONTROLLER) DEFINED BY 
GFC BASED ON QSAE MODELS 

• INCORPORATED THE CONTROLLER IN DYNAMIC EQUATIONS OF 
MOTION IN TO THE BOEING (SEATTLE) P-K FLUTTER ANALYSIS 
PROCESS 

• ESTABLISHED OPEN LOOP AND CLOSED LOOP FLUTTER SOLUTIONS 
FOR THE PITCH CONTROLLER 

• ESTABLISHED WITH FLIGHT CONTROLS TO VERIFY CONCURRENCE 
BETWEEN P-K FLUTTER SOLUTION RESULTS & STATE-SPACE MODEL 
RESULTS 
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DYNAMICS MODEL 


• DITS TCAY MODEL OPTIMIZED FOR STRENGTH AND FLUTTER 

• MASS CONDITIONS CONSIDERED - Ml, M02, ML2, MCF, MCI, MT1, MT4 

• UNSTEADY AERODYNAMICS 

DOUBLET LATTICE TYPE METHOD FOR M = .24, .40, .65, .80, .90, .95 
A502 AERO FOR M = 1.20, 2.0, 2.6 

REDUCED FREQUENCIES ‘k’ FROM 0.0 TO .0127 ( 21 VALUES ) 

• COMMON MODEL FOR ALL OF THE FOLLOWING 

DESIGN INTEGRATION TRADE STUDIES 
AEROELASTIC DESIGN STUDIES 
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TCA Finite Element Model (FEM) 
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TCA Elfini Flutter Aerodynamic Meshes 
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TCA Elfini Flutter Computational Grid 


Shaded regions denote zones within which the Elfini basis deformation 
shapes are defined 

Coincident zones are defined on each aerodynamic mesh 
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Generalized airforces are computed in terms of basis shapes 
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High Speed Civil Transport 


FLUTTER EQUATIONS OF MOTION (APEX P-K PROGRAMS : 
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PITCH COMTROL LAWS (SIMPLIFIED GAMMA DOT- V CONTROLLERS : 
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SIMPLIFIED ACTUATOR MODEL: 
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TCA CONFIGURATION 
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FLIGHT ENVELOPE 
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Longitudinal Control Law Overview 


HSCT Flight Controls and the ASE Working Group 


Longitudinal Control Law Overview: jV 


• Present Control Laws 


- The throttle controls speed. 

- The elevator controls flight path. 



• j/V Control Law 

- The throttle controls the total energy. 

- The elevator controls how the energy is split between flight path and speed. 
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Longitudinal Control Law Overview: 7 V (Cont.) 


• Basic concept of j/V is to control total energy 

• Total energy is sum of potential and kinetic 

1 2 

et = mgh -| — mV 

2 

• Specific energy rate (normalized by weight and velocity): 

T — D V 

E t = = 7 H 

mg g 


(1) 

( 2 ) 




f~ I <£f U A £ (3 2. Lo rJ 6; /r U i) I i do «JT A. oL 1—4 t*J O £■ 4 V/ / <£ UJ 



Ccc Afty 



HSCT Flight Controls and the ASE Working Group 


Longitudinal Control Law Overview: 7 V (Cont.) 


• Along flight path, required trim thrust is: 

^~req = mgE t + D (3) 

• Assuming drag variation is slow, T req is proportional to specific energy rate (i.e. 
throttles control rate of energy addition) 

• Use throttle feedback law: 

^ A (V XP + E te = (k t p + (4) 
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Longitudinal Control Law Overview: 7 V (Cont.) 


• Forces resulting from elevator movement are essentially conservative — > no 
in total energy. 

• Elevators do not significantly change aircraft’s total energy state. 

• Elevators serve mainly to redistribute energy between potential and kinetic. 

• Use elevator to control energy distribution error = y — 7. 

• Elevator feedback law: 

S E = (ksp + E de = (i K SP + - 7^ 
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Simplified jV Control Law Block Diagram 
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Simplified 7V Control Law Equations 
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Total Closed Loop System 
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State-Space Modeling Methodology 
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FC DASE Modeling Methodology 



Flutter Equation 


Ms 2 + Cs + K 

- qA ua (s)] X(a) = $ T GU(s) + q [YgWAgy 

(s)+'Z g (a)A gz (S)] 

(10) 

where 

A gy (s) 4 Y,A igy {s)e-^ s 

i= 1 

(11) 
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FC DASE Modeling Methodology 


• Flight Controls derives state space models from Structures model (File 41) 

- Automated routine to convert Elfini File 41 to Matlab format. 

- Collection of Matlab m-files used to: 

1. Fit unsteady aero using Roger RFA 

2. Introduce structural viscous damping 

3. Define control inputs 

4. Produce sensor outputs 

5. Return state space A, B, C, D matrices 

• Matlab routines also developed to perform: 

1. Frequency domain response analysis 

2. Flutter analysis {V — g diagrams) 
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FC DASE Modeling Methodology 


Unsteady Aero Approximations 
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FC DASE Modeling Methodology 

State Space Equations of Motion 

•^ss — -^ss^ss “E B ss u ~J~ B SS gW (16) 

X ss = ® X r l r 2 r Uj3 T'gy 1 r gy 2 • * * r gyn j3 r gzi r gz 2 • • • r gz r . 

( 17 ) 

W T = [ jjg Vg Zg Zg ] (18) 
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FC DASE Modeling Methodology 


Torque Inputs vs Prescribed Motion 

• Equation (16) uses generalized internal/external forces as driving inputs. 

— Actuator force model required to drive system. 

— Accounts for full two-way coupling of surface-into-vehicle and vehicle-into-surface 
dynamic coupling. 

— Surface motion is an output. 

• Prescribed surface motion commonly used when actuator models not well defined. 

— Actuator displacement model required to drive system. 

— Accounts only for surface-into-vehicle dynamic cou pi ing. 

— Equation (16) may be transformed to prescribed motion, making required hinge 
moment an output. 
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Flight Controls 


State-Space Model Generation Process 
' Flight Controls 

1. Convert Elfini File 41 data to Matlab format 

2. Define State-space Model Build Specification 

• Model file name 

• Added viscous structural damping 

• Flight condition (Mach and altitude) 

• Retained modes 

• Unsteady aero fitting parameters 

- Number of lag states 

- Force fit at zero frequency 

— Use inverse frequency weighting 

- Omit apparent mass term 

• Define sensor and actuator support node numbers 

3. Determine true airspeed and density 

4. Form basic A, B, C, D representation 

• Calculate viscous damping matrix 

• Determine indices of retained modes 

• Strip-out unretained modes 

• Calculate dynamic pressure 

• Determine frequencies of lag filters 

— Frequencies are evenly spaced between minimum and maximum fc-values 

• Form state-space A matrix 

- Fit unsteady aero using linear weighted least squares 

- Combine aero fitting coefficients with mass, stiffness, and damping matrices to form 
state-space A matrix 

— Correct A matrix for accelerations due to plunge position (should be zero). 

— Correct A matrix such that accelerations due to rigid body plunge velocity agree 
with accelerations due to rigid body pitch position, the latter being the truth model. 

• Form state-space B matrix 

— Collect mode shapes for assumed modes, or actuator degress of freedom (Identity 
matrix) 

- Collect mode shapes for actuator support structure 

- Use finite differencing, if necessary, to obtain modal slopes 

- Calculate difference between mode shapes for the actuator DOFS and the actuator 
support structure 

- Form state-space B matrix using the mass matrix, apparent aerodynamic mass, and 

the mode shape differences. 
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Flight Controls 


• Form C matrix ( 

— Include rigid-body mean-axis states 

- Include nodal displacements 

— Include assumed mode deflections 
— Strip out unretained modes 

- Calculate C matrix for y 
— Calculate C matrix for y 

— Form total C matrix for y, y, and y 

• Form D matrix 

- Include direct feed-through elements from y 

• Convert A, B, C, D representation to sparse arrays 

5. Add additional sensors to C and D matrices 

• Pitch and pitch rate sensors at various body stations 

• Load factor sensors at various body stations 

• Velocity sensors at various body stations 

• Assemble new combined C and D matrices 

6. Transform coordinates and convert units 

• Define coordinate transformations 

• Define unit conversions 

• Form transformations matrices for state vector, measurement vector, and control vector. 

• Perform transformation 

7. If desired, convert from torque inputs to prescribed motion inputs 

8. Done 
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PART II - Boeing Long Beach Report 


Introduction 

A closed-loop ASE flutter analysis of the TCA configuration with symmetric boundary 
condition has been performed using MSC/NASTRAN. While the primary goal of this 
study was to validate the MSC/NASTRAN ASE analysis process and not to generate ASE 
stability results for the TCA configuration, several interesting results were nonetheless 
generated. 

MSC/NASTRAN ASE Analysis Process 

MSC/NASTRAN is an attractive package for ASE analysis since the control laws can be 
directly included with the structural and aerodynamic models in a single ASE model, 
without requiring the intermediate steps of generating vibration modes, generalized AIC’s, 
and the like. In a DITS-type sizing process, the transfer functions could conceivably be 
included directly in the flutter screening analyses as well as in the flutter constraint 
calculation in the strength/loads/flutter sizing process. 

While MSC/NASTRAN allows the inclusion of control laws directly in the flutter analysis, 
there are several things that must be kept in mind, and issues that must be overcome. 

1. The structural FEM must have rigid-body (i.e. zero frequency) 
vibration modes for each of the control surfaces to be used in the 
analysis. For this study, the elevator was the only surface used, so the 
model had four rigid body modes (pitch, plunge, fore-aft, and elevator 
pitch). 

2. The FEM must have degrees of freedom that directly correspond to the 
sensors of the ASE model. For example, if an acceleration at the pilot 
station is a desired input, then a structural node must exist at the pilot 
station. 

3. A structural node must also be created for each actuator, and MPC’s 
defined relating the actuator node deflection to the physical deflection 
of the associated control surface 

4. MSC/NASTRAN has no built-in gain scheduling capability. If gain 
scheduling is required, then a separate transfer function must be defined 
for each speed point, and each speed point must be analyzed in a 
separate subcase. 

Addressing these issues required several minor modifications to the finite element model 
as well as writing a simple preprocessor code to automate the gain scheduling task. 



Closed Loop Flutter Analysis Results 


In validating the closed loop MSC/NASTRAN flutter analysis two mass conditions, MO-2 
and MT-1, were investigated. The open loop flutter results for some of the critical modes 
of these two mass conditions are presented in Figures 1 and 2. It can be seen that the 
open-loop aircraft is aeroelastically stable within the flight envelope. 

The Gamma dot V control was included in the ASE model for this validation process. 
That control law, including actuator dynamics, takes the form: 
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^ ELEV ~~ 25 STAB 


The Gamma dot V control law drives both the elevator and the horizontal stabilizer. In the 
analysis presented in this memo the stabilizer deflection was not included, and the control 
law only drives the elevator. Since the elevator gain is twice that of the horizontal 
stabilizer, it is believed this modification would not affect the conclusions of this analysis. 
The sensor configuration used in these analyses were located near the center of mass of 
the aircraft. 


The closed loop flutter analysis results for mass conditions MO-2 and MT-1 are presented 
in figures 3 and 4 respectively. It can be seen that both conditions exhibit damping in the 
unacceptable range. 

Raj Nagaraja presented closed loop flutter results for the MO-2 condition at the HSR 
Airframe Technical Review in February 1998. The results presented in this memo 
correlate well with those results. This correlation validates the MSC/NASTRAN closed- 
loop flutter analysis. 


Possible Modifications to the “Gamma dot V” Control Law 


A further investigation was performed to try and devise a filtering scheme that would 
reduce the effects of the control law on the flexible aircraft. It was determined that the 
main cause of the instabilities was the pitch rate feedback (q), and therefore, only this loop 
was filtered. The control employed in the analysis then became: 
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Where NTF represents the filter transfer function. Initially a single notch filter was 
employed at the instability frequency. It was found that such a filter would reduce the 
initial instability, but would cause sufficient phase shift to drive another mode unstable. 
Subsequently a double notch configuration was devised that addressed all aeroelastic 
instabilities. These filters consist of 14db notch filters at 1.8Hz and 3.6Hz and combined 
to form the transfer function: 
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S 2 +3.1645 + 127.69 
S 2 +15.825 + 127.69 


^ S 2 +6.3285 + 5 10.76 ^ 
S 2 +31.645 + 510.76 


) 


A bode plot of this filtering scheme is presented in figure 5. Aeroservoelastic analysis 
results incorporating this filtering scheme are presented in figures 6 and 7 for mass 
conditions MO-2 and MT-1. It can be seen in Figure 7 that at the low speed end of the 
flight envelope a humped mode crosses over into the region of unacceptable damping. 
Further refinements of the notch filters would probably eliminate this instability. 

This filtering scheme and the accompanying analyses are presented as a possible solution 
to the aeroservoelastic instability problem. It is not meant to be the final or only solution. 
Some further filter refinement and analysis is still required to: 

1. Completely eliminate all aeroelastic instabilities within the flight envelope, for 
all mass and flight conditions 

2. Determine the filters’ effects on the aircraft’s flying qualities. 

















Frequency vs Velocity 
.95M Mass Condition MT1 
Baseline Gamma-dot V Control Law 



Damping vs Velocity 
.95M Mass Condition MT1 
Baseline Gamma-dot V Control Law 



Figure 4: Closed Loop Flutter (Baseline Gamma-dot- V) Analysis Results, Mass 

Case MT-1 at Mach 0.95. 
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Figure 5: Bode Plot of Double-Notch Filter Used to Eliminate Aeroservoelastic 

Instabilities. 
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Frequency vs Velocity 
.95M Mass Condition M02 
Gamma-dot V Control Law with Notch Filters 
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Damping vs Velocity 
.95M Mass Condition M02 
Gamma-dot V Control Law with Notch Filters 
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Figure 6: Closed Loop Flutter (Ganima-dot-V with Notch Filter) Analysis Results, 

Mass Case MO-2 at Mach 0.95. 





Figure 7: Closed Loop Flutter (Gamma-dot-V with Notch Filter) Analysis Results, 

Mass Case MT-1 at Mach 0.95. 
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